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tum. This mechanism will provide a clear experimental signal for Higgs production 
due to the small background in this kinematic region. The key assumption under- 
lying our analysis is the presence of intrinsic heavy flavor components of the proton 
bound state, whose existence at high light-cone momentum fraction x has growing 
experimental and theoretical support. We also discuss the implications of this picture 
for exclusive diffractive quarkonium and other channels. 
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1 Introduction 



A central goal of the Large Hadron Collider (LHC) being built at CERN is the 
discovery of the Higgs boson, a key component of the Standard Model, and whose 
discovery would constitute the first observation of an elementary scalar field. A 
number of theoretical analyses suggest the existence of a light Higgs boson with a 
mass Mh < 130 GeV. 

In this paper we propose a novel mechanism for hadronic Higgs production, in 
which the Higgs is produced with a significant fraction of the projectile momentum. 
The key assumption underlying our analysis is the presence of intrinsic charm (IC) 
and intrinsic bottom (IB) fluctuations in the proton bound state [1, 2], whose exis- 
tence at high X as large as x ~ 0.4 has a substantial and growing experimental and 
theoretical support. Clearly, this phenomenon can be extended to the consideration of 
intrinsic top (IT). A recent review of the theory and experimental constraints on the 
charm quark distribution c{x, Q'^) and its consequences for open and hidden charm 
production has been given by Pumplin [3]. The presence of high x intrinsic heavy 
quark components in the proton's structure function will lead to Higgs production 
at high xp through subprocesses such as gh ^ Hb; such reactions could be particu- 
larly important in MSSM models in which the Higgs has enhanced couphngs to the b 
quark [4]. 

The virtual Fock state \uudQQ > of a proton has a long lifetime at high energies 
and can be materialized in a collision by the exchange of gluons. The heavy quark 
and antiquark can then coalesce to produce the Higgs boson at large xp — Xc-\- Xc- 
This Higgs production process can be inclusive as in pp — > HX, semi-diffractive 
pp — > HpX, where one of the projectile protons remains intact, or exclusive diffractive 
PP ~^ pHp, where the Higgs can be reconstructed from the missing mass distribution. 
In each case the Higgs distribution can extend to momentum fractions xp as large as 
0.8, reflecting the combined momentum fractions of the heavy intrinsic quarks. 

Perhaps the most novel production process for the Higgs is the exclusive diffractive 
reaction, pp — >■ p + H + p [5], where the + sign stands for a large rapidity gap 
(LRG) between the produced particles. If both protons are detected, the mass and 
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momentum distribution of the Higgs can be determined. The TOTEM detector [6] 
proposed for the LHC will have the capability to detect exclusive diffractive channels. 
The detection of the Higgs via the exclusive diffractive process pp ^ p + H + p, has 
the advantage that it does not depend on a specific decay mechanism for the Higgs. 
The branching ratios for the decay modes of the Higgs can then be individually 
determined by combining the measurement of cr{pp p + H + p) with the rate for 
a specific diffractive final state Bf a{pp p + H^f + p). This is in contrast to the 
standard inclusive measurement, where one can only determine the product of the 
cross section and branching ratios Bf a(pp — > H^fX). 

The existing theoretical estimates for diffractive Higgs production are based on 
the gluon-gluon fusion subprocess, where two hard gluons couple to the Higgs {gg 
H) [5]. A third gluon is also exchanged in order that both projectiles remain color 
singlets. Pcrturbativc QCD then predicts a[pp — >■ p-\-H +p) ~ 3 fb for the production 
of a Higgs boson of mass 120 GeV at LHC energies, with a factor of 2 uncertainty 
[5] . Since the annihilating gluons each carry a small fraction of the momentum of the 
proton, the Higgs is primarily produced in the central rapidity region. 

In this paper we will specifically consider the exclusive diffractive production re- 
action pp — > p-\-M-\-p, depicted in Fig. 1, where M stands for J/il>, Xc Xb, or H. 
The final state M will be produced in the projectile proton's fragmentation region 
with a significant fraction of the incident proton's momentum, since the sum of the 
momenta of two heavy quarks contribute to the momentum of M. This has an im- 
portant advantage of providing a distinctive signal with relatively small background. 
This production process is analogous to the positron-antiproton coalescence reaction 
by which anti- hydrogen was first detected [7, 8, 9]. 



2 Intrinsic Heavy Quarks 

n; Xi, k±^, \i ) of the QCD Light- 



The proton eigenstate \p) = ^n=3'^ri{xi, k±^, Xi) 
front Hamiltonian H^p^ can be expanded at fixed light-front time r = t + z/c 
as a superposition of quark and gluon Fock eigenstates of the free Hamiltonian 
^LF '■ \p) — \uud), \uudg),---, including, in particular, a "hidden charm" Fock 



3 




Figure 1: The exclusive difFractive production of J/ip, Xc T", Xb, Z'^ oi H, the Standard 
Model Higgs. 

component | uudcc) . The fact that the hadronic eigenstate has fluctuations with 
an arbitrary number of constituents is a consequence of quantum mechanics and 
relativity. The | uudcc) Fock state arises in QCD not only from gluon splitting 
which is included in DGLAP evolution, but also from diagrams in which the heavy 
quark pair is multi-connected to the valence constituents. The latter components 
are called "intrinsic charm" (IC) Fock components. The frame-independent light- 

— * 

front wave functions ipnixi, k±^, Aj) describe the constituents in the hadron with mo- 
menta — XiP'^,p±i — XiP± + k±i and spin projection Aj. Here Yl^^iXi — 1, and 
^^^^ = 0. The hght-cone momentum fractions Xi = p+/P+ = {pf + Pf) / + P^) 
are boost invariant [10]. 

It was originally suggested in Refs. [1, 2] that there is a ~ 1% probability of IC 
Fock states in the nucleon; more recently, the operator product expansion has been 
used to show that the probability for Fock states in hght hadron to have an extra 
heavy quark pair of mass Mq decreases only as ^qcd/^q non-Abelian gauge the- 
ory [11]. In contrast, in the case of Abelian QED, the probability of an intrinsic heavy 
lepton pair in a light-atom such as positronium is suppressed by fif^f^^/M^, where fibohr 
is the Bohr momentum. The quadratic QED scaling corresponds to the dimension-8 
Euler-Heisenberg effective Hamiltonian F'^/Mf for light-by light scattering mediated 
by heavy leptons. Here F^j^ is the electromagnetic field strength. In contrast, the cor- 
responding effective Hamiltonian in QCD /Mq has dimension 6. This difference in 



power behavior provides a remarkable discriminant between non-Abelian and Abelian 
theory. 

The maximal probability for an intrinsic heavy quark Fock state occurs for mini- 
mal off-shellness; i.e., at minimum invariant mass squared Ai^ = Y17=3i''^l + ^iJ/^i- 
Thus the dominant Fock state configuration is Xi oc m±i where m^_|_ = + /c^^^; i.e., 
at equal rapidity. Since all of the quarks tend to travel coherently at same rapidity 
in the | uudQQ) intrinsic heavy quark Fock state, the heaviest constituents carry the 
largest momentum fraction [1, 2]. Models for the intrinsic heavy quark distributions 
ci{x, Q"^) and hi{x, Q^) predict a peak at x ~ 0.4. Thus the intrinsic heavy quarks are 
highly efficient carriers of the projectile momentum. 

Intrinsic charm also leads to new, competitive decay mechanisms for B decays 
which are nominally CKM-suppressed [12] and in explaining the J/ip pn puz- 
zle [13]. Furthermore, it has been found that intrinsic bottom could even contribute 
significantly to exotic processes such as neutrino-less ii~ — e~ conversion in nuclei [14]. 

3 Relevant Experimental Facts 

The most direct test of intrinsic charm is the measurement of the charm quark dis- 
tribution c(x, Q"^) in deep inelastic lepton-proton scattering ip i'cX. The only 
experiment which has looked for a charm signal in the large Xhj domain is the Euro- 
pean Muon Collaboration (EMC) experiment [15], which used prompt muon decay in 
deep inelastic muon-proton scattering to tag the produced charm quark. The EMC 
data show a distinct excess of events in the charm quark distribution at x^j > 0.3, at 
a rate at least an order of magnitude beyond predictions based on gluon splitting and 
DGLAP evolution. Next-to-leading order (NLO) analyses [16] show that an intrinsic 
charm component, with probabihty of order 1%, is needed to fit the EMC data in 
the large Xbj region. This value is consistent with an estimate based on the operator 
product expansion [11]. Clearly it would be very valuable to have additional direct 
measurements of the charm and bottom structure functions at large x. 

An immediate consequence of intrinsic charm is the production of charmonium 
states at high xp — Xc + Xc in high energy hadronic collisions such as pp ^ J/ipX. 
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The c and c in the IC Fock state | uudcc) can be materiahzed by gluon exchange 
as a color-singlet pair which coalesces to a high xp low pr quarkonium state. The 
internal color structure of the Fock state is important. The effective operator in the 
non-Abelian theory predicts that the charm quark pair is dominantly a color-octet 
8c- The color octet {cc)sc is then converted to a high x color-singlet (cc)i^ state via 
gluon exchange with the target; it then couples to the color-singlet quarkonium state. 
Note that the J/ip can be produced this way only from the component of IC which 
is symmetric, relative to a simultaneous permutation of spatial and spin variables. 

Comprehensive measurements of the pA — > J/ t/jX and nA — > J/ t/jX cross sections 
have been performed by fixed target experiments, NA3 at CERN [17] and E886 at 
FNAL [18]. According to the arguments in Refs. [19, 20, 21], the IC contribution is 
strongly shadowed, thus accounting for the observed nuclear dependence of the high 
xp component of the J/ip hadroproduction. It is also important to consider effects 
coming from energy conservation. Multiple interactions in a nucleus can resolve higher 
Fock components of the projectile hadron compared to interaction with a free proton 
target. Therefore, energy sharing between the projectile partons imposes more severe 
restrictions on production of energetic particles leading to nuclear suppression at large 
xf [22]." 

The materialization of the intrinsic charm Fock state also leads to the production 
of open-charm states such as A{cud) and D^(cd) at large xp- This may occur either 
through the coalescence of the valence and charm quarks which are co-moving with 
the same rapidity, thus producing a leading particle effect or via hadronization of 
the produced c and c. As shown in Refs. [23, 24], a model based on intrinsic charm 
naturally accounts for the production of leading charm hadrons in pp — > DX and 
pp — > AgX as observed at the ISR [25] and also at Fermilab [26, 27]. We also note 
that it is also possible to construct Regge models which give similar xp behavior as 
the IC approach. 

The diffractive cross section a{pp AcX +p), at ^/s = 63 GeV was measured at 
the ISR to be is of order 10 to 60 /x6 [28]. This result seems to be in contradiction with 
findings of Fermilab experiments searching for diffractive charm production [29, 30]. 
The E690 experiment [29] observed the diffractive channel a{pp D*pX) ~ 0.2 jib 
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at ^/s — 40 GeV. Their results lead to the diffractive charm production cross section 
^difficcy^ — 0.6 — 0.7 /lb which is about two order of magnitude smaller than the 
cross section measured at ISR [28] . This result agrees with the upper hmit found for 
coherent diffractive production of charm, pSi ccX+Si, in the E653 experiment [30] 
at Fermilab. However, forward charm production is most likely strongly suppressed 
in a nuclear target as is the case for light hadrons. If one extrapolates to pp collisions 
assuming an A^/^ dependence [31], the upper hmit is Udiff^cc)'^ < 1 jih. The ISR 
signals for forward charm production are thus not necessarily inconsistent with the 
fixed target experiments considering the large differences in the available center of 
mass energy, as well as the nuclear target suppression. 

It should be noted that diffractive charm production via elastic scattering of 
the projectile plus gluon radiation also leads to the right order of magnitude of the 
cross section. Indeed, the cross section for diffractive gluon radiation (via the triple- 
pomeron term) in pp collisions is known from data, a^f ~ 4 mb. The production of 
a cc pair via a radiated gluon brings extra factors from the coupling ct^ 0.2 and 
from the gluon propagator (m^/Mcc)^ where rUg ~ 0.8 GeV is the effective gluon 
mass [32]. Thus, one obtains an estimate for the singly diffractive charm production 
cross section a{pp — > ccpX) 4 in good agreement with the magnitude of the 
data, at least in the central rapidity region. However, the shape of the empirical Ac 
distribution at large xp \s not readily accounted for this model. 

One can produce the K{hud) at high xp in inclusive pp collisions, through the 
materialization of the intrinsic bottom Fock state | {uud)^^^,{bb)g^y The cross section 
for forward open bottom production relative to open charm is reduced by the relative 
IC/IB probability factor ml/ ml ~ 1/10. Evidence for the forward production of the 
A;, in pp — > AfjX at the ISR was reported in Ref . [33] . 

The existence of rare double-IC Fock state fluctuations in the proton, such as 
\uudcccc > can lead to the production of two J/ip's [34] or a double-charm baryon 
state at large xp and small Pt- Double- J/-?/' events at a high combined xp > 0.5 
were in fact observed by NA3 [35]. The observation of the doubly-charmed baryon 
S+(3520) with mean < xp >~ 0.33 has been reported by the SELEX collaboration 
at FNAL [36]; the presence of two charm quarks at large xp has, indeed, a natural 
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IC interpretation. 



4 Intrinsic Heavy Quarks and Exclusive Diffrac- 
tive Production 

We now investigate the implications of IC and IB for exclusive diffractive production 
processes pp — > pMp at large xp where M is a charmonium state, boson, or the 
Higgs. We first explain, using Fig. 2, how the exclusive diffractive channels shown in 
Fig. 1 arise with the required color structure in the final state. As noted above, we 
shall assume that the projectile (upper) proton has an approximate 1% probability 
to fluctuate to an IC Fock component with the color structure |[Mwci]8c[cc]8c >• This 
virtual state has a long coherence length in a high energy collision oc s/M.'^Mp, 
where Ai is the total invariant final-state mass. In a pp collision, two soft gluons 
must be exchanged in order to keep both protons intact and to create a rapidity 
gap, mimicking pomeron exchange. The two gluons couple the target nucleon to the 
large color dipole moment of the projectile IC Fock state. For example, as shown in 
Fig. 2, one of the exchanged gluon can be attached to the d valence quark spectator 
in |[MM(i]8(j[cc]8c >, changing its color, and the other one can be attached to the c, 
also changing its color. The net effect of this color rearrangement is the same as 
single-gluon exchange between the two color-octet clusters. The cc and the uud can 
thus emerge as color singlets because of the gluonic exchange. The [cc] can couple 
to the J/'0, or to a Z° or to a H. Meanwhile the color-singlet uud gives rise to the 
scattered proton, thus producing the two required rapidity gaps in the final state. 
Notice that the xp distribution of the produced particle is approximately the same 
as the distribution of the [cc] inside the proton. As we shall discuss below, the sum of 
couplings of the gluon to all of the quarks, as dictated by gauge invariance, brings in 
a form factor which vanishes at zero momentum transfer, thus giving an important 
suppression factor. 
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4.1 The cross section 

The cross section of exclusive diffractive production of the Higgs, pp Hp + p, can 
be estimated in the hght cone (LC) dipole approach [37]. The Born graph for this 
process is shown in Fig. 2. As discussed above, we shall assume the presence in the 



H^) 



P = 
mr,P;z) 




Figure 2: The two-gluon exchange diagram for the Higgs exclusive pro- 
duction 



proton of an intrinsic charm (IC) component, a cc pair, which is predominantly in 
a color-octet state, and which has either a nonperturbative or perturbative origin. 
In the former case this heavy component can interact strongly with the 3q valence 
quark component. Such nonperturbative reinteractions of the intrinsic sea quarks in 
the proton wavefunction can lead to a Q{x) ^ Q{x) asymmetry as in the AK model 
for the ss asymmetry [38, 39]. As in charmonium, the mean cc separation should be 
considerably larger than the transverse size l/rric of perturbative cc fluctuations. For 
instance, if the binding potential has the oscillator form, the mean distance is 



where uj ~ 300 MeV is the oscillation frequency. Alternatively, the IC component 
can be considered to derive perturbatively from the minimal gluonic couplings of the 
heavy quark pair to two valence quarks of the proton; this is likely the dominant 
mechanism at the largest values of Xc [21]. In this case the transverse separation of 
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the cc is controlled by the energy denominator, (r^,) = i/iril, and is much smaller 
than the estimate given by Eq. (1). 

In accordance with the notation in Fig. 2 the two protons in the CM frame are 
detected with Feynman momentum fractions Xi and X2 and transverse momenta pi 
and P2 respectively. Correspondingly, the produced Higgs carries longitudinal mo- 
mentum (^2 — Xi) and transverse momentum pn — —{pi +^2)- We assume the Higgs 
to be heavy, over 100 GeV, then xi and X2 turn out to be tightly correlated in this 
reaction. Indeed, the effective mass squared of the H — pi pair reads, 

1-Xi Xi ^ ^ ^ 

(We assume the equivalence of Feynman x and the corresponding fractions of the 
light-cone momenta, which is an accurate approximation at large x.) Due to the 
form factors of the two protons, neither transverse momenta, pi^2, can be much larger 
than a few hundred MeV and therefore they, together with ph and the proton mass, 
can be safely neglected in Eq. (2). This could be incorrect at very small values of 
Xi ~ rrip/M'^, but we will show that the Xi distribution sharply peaks at Xi ~ 0.25. 
Then, employing the standard relation M'^ / s = I—X2, we arrive at the simple relation, 

(1-Xi)(l-X2) = ^ . (3) 

s 

The diffractive cross section has the form, 

da{pp-^ppH) 1 - -m2 rA\ 

dx2d^p,d-p2 = (l-a;2)16vr^ '^^^^'^^'^^^' ' 

where the diffractive amplitude in Born approximation reads, 

8 [ (Pq d^k ^ ^ 

A{x2,Pi,P2) = -^-j^ j d^Q — —as{q^)as{k'^)5{q + p2 + k)5{k-pi-Q) 

X (1 - e-'^-^ $J(p)e^^-^72 (^i _ e-ife-p-^ ^^(^^ ^) ^iQ-R (5) 
Here r, p, z) is the light-cone wave function of the IC component of the projectile 
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proton with transverse separations R between the cc and 2>q clusters, r between the c 
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and c, Q is the relative transverse momentum of the ?>q and cc clusters in the projectile 
and p is the transverse separation of the quark and diquark which couple to the final- 
state proton p2- The density |$p(r)p is the wave function of the target proton which 
we also treat as a color dipole quark-diquark with transverse separation r. (The 
extension to three quarks is straightforward [37]). The fraction of the projectile 
proton light-cone momentum carried by the cc, 2; ~ 1 — xi. This wave function is 
normalized as, 

1 

' = Pic. (6) 


where Pjc is the weight of the IC component of the proton, which is suppressed as 
and is assumed to be Pic ~ 1%- The amplitudes H{f) and ^p{p) denote 
the wave functions of the produced Higgs and the outgoing proton, respectively, in 
accordance with Fig. 2. 

The phase factors in Eq. (5) correspond to different attachments of the exchange 
gluons to quarks in Fig. 2. Thus, attaching the gluon either to the c, or to the c quarks 
one gets factor [exp(ig • r/2) — exp(— • r/2)]. An analogous factor corresponding to 
the second gluon coupling to the proton pi is also included in Eq. (5). The transverse 
coordinates of the quark and diquark in the target proton are t/2 and —r/2 (relative 
to its center of gravity). The phase factor in the square brackets in Ref. (5) thus 
includes two terms corresponding to attachment of the exchanged gluons to the same 
or different valence quark or diquark in p2. 

In order to advance the calculations further, we will take the following steps: First, 
we assume a factorized form of the proton wave function, 

^>,{R, f, p, z) = ^>ic{R. z) *3c(r-) *3,(p) . (7) 

Here and '^^q are the cc and 3g wave functions normalized to unity, whereas 
ic{Ri z) is the wave function describing the relative motion of the cc and 3g clusters, 
where z is the fraction of the longitudinal momentum carried by the cc. This wave 
function is normalized as. 
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dz I (fRcfrcfp 



^p{R,f,p,z) 



where Pic{z) is the ^^-distribution of cc, related to the xi distribution of the produced 
protons, since with a very high precision z — 1 — Xi — Mjj/ s{l — X2) (unless Xi is as 
small as xi ~ 2mp/ . 

We will perform the calculations in Eq. (5) only for forward diffraction, i.e. p2 = 0, 
q — —k, and we assume for the Pomeron the typical Gaussian t-dependence {t = — 



da 



so the t-integrated cross section then reads, 

da TT da 



oc e-^(^'>^ , (9) 



(10) 

P2=0 



d'^Pidx2 B{s') d'^pi d'^p2 dx2 

Here the slope B{s') ~ So + 2q;^ ln{s'/M^), where = 4 GeV~^ a'p = 0.25 GeV~^ 
s'/M^ = s/M^ and Mq = 1 GeV. 

The next step is to replace the two-gluon proton vertex, represented by the integral 
over -fin Eq. (5), by the unintegrated gluon density, J^{x,k!^) = dG{x, k'^) /d{\nk'^), 
where G{x, Q"^) = x g{x, Q"^). This preserves the infra-red stability of the cross section, 
since JF vanishes at A;^ — * 0. The phenomenological gluon density fitted to data 
includes by default all higher order corrections and supplies the cross section with an 
energy dependence important for extrapolation to very high energies. One can relate 
the unintegrated gluon distribution to the phenomenological dipole cross section fitted 
to data for ^2(2:, Q"^) from HERA, as was done in Ref. [40], 



Hx, k') = k' Rlix) exp \-\ Rlix) k' 
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The problem is the extrapolation to the small virtualities A;^ typical for the process 
under consideration. The Bjorken variable is not a proper variable for soft reactions; 

therefore we use the parametrization from Ref. [32] adjusted to data for soft interac- 
tions. Then Ro{x) in Eq. (11) should be replaced by Rq{s') = 0.88 fm x (s'/sq)"'^^^ 
with A = 0.28, So = 1000 GeV^ and ao ^ ao{s') = aZi-s') [1 + 37?g(s')/8(r2^)^], where 
(7^^(s') — 23.6 mb x (s'/so)° °^ is the Pomeron part of the 7rp total cross section. The 
energy variable s' is related to the rapidity gap between the two protons in the final 
state, controlled by X2, 
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Finally, combining all the above modifications and performing the pi-integration 
in Eq. (5), we arrive at, 



da^^{pp ppH) 327r Pic{z) 



dx2 9B{s'){l-X2) 

X 



X J d'p '^l{p)e-''-^/' (l - e*'^) ^3,(p) 



(13) 



Here 



z^[x^ + ^(xf)2 + 4Mys]/2 ^xf^l-xi^ Ml/s{l - X2) . (14) 

This relation receives sizeable corrections only at very small Higgs Feynman Xp ~ 
2Mh/\/s- Notice that the expansion of the exponentials in Eq. (13) contains only 
odd powers of k ■ r and k ■ p. This signals a change of orbital momentum of the quark 
configurations participating in the onc-gluon exchange process. In order to obtain a 
nonzero result of the integration over r, either the initial or the final cc wave function 
must contain a factor V^, i.e. it must be a P-wave. Since we assume that the Higgs 
is a scalar, its cc component must be in a P-wave state, while the primordial cc in 
the projectile IC state should be in an S-w&ve. This is vice versa for the proton pi. 
the final \3q) system is in an 5"- wave, but ^zq{p) must be a P-wave. 

Notice that both the scalar Higgs and x states may be produced from the same 
IC component of the proton containing S'-wave cc. However, the production of J/ip., 
T, requires an IC component containing a P-wave cc, which is presumably more 
suppressed. 

The P-wave LC wave function of Higgs in impact parameter representation is 
given by the Fourier transform of its Breit-Wigner propagator: 

^ (15) 



H[r) = I — rric X cr X - 

2n r 



IT 

eY,{er)--THMHYo{er) 



Here is the Fermi constant, x X spinors for c and c respectively and 

e2 ^ a{l - a)Ml - ml , (16) 

where a is the fraction of the LC momentum of the Higgs carried by the c-quark. The 
functions 1^0,1(2^) iii Eq. (15) are the second order Bessel functions and Vh is the total 
width of the Higgs. Assuming Tq <^ Mh, we neglect the second term in Eq. (15). 
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The LC wave function Eq. (16) assumes that the Higgs mass is much larger than 
the quark masses, which is probably true for charm and bottom. However, it is quite 
probable that for top-antitop in the Higgs 2mt > Mh, then the wave function is 
different, 

H-ttir) = ^^"^^ mtX^X-^tKi{etr) , (17) 

ZTT r 

where Ki{x) is the modified Bessel function and 

e^t ^ml-a(l-a)Ml . (18) 

The probabilities computed from the wave functions Eqs. (15) and (18) require 
regularization in the ultraviolet limit [41, 42], as is the case of the qq wave function of a 
transverse photon. Such wave functions are not solutions of the Schrodinger equation, 
but are distribution functions for perturbative fluctuations. They are overwhelmed by 
very heavy fluctuations with large intrinsic transverse momenta, or vanishing trans- 
verse separations. Such point-like fluctuations lead to a divergent normalization, but 
they do not interact with external color fields, i.e. they are not observable. All the 
expressions for any measurable quantity, including the cross section, are finite. 

As we have discussed, the IC wave function can be modeled as a nonperturbative 
5 — quark stationary state |3gcc), or as a perturbative fluctuation \3q) \3qcc). 
Correspondingly, the cc wave function within the Fock state will be assumed to be a 
linear combination of nonperturbative and perturbative distribution amplitudes, 

*.c(r) = P *r (^1 + a/i^ "^ti^ ■ (19) 

The parameter f3 which controls the relation between the nonperturbative and pertur- 
bative IC contributions, is such that < f3 < 1. The nonperturbative wave function 
should be an S-wave solution of the Schrodinger equation. Assuming oscillator po- 
tential we get, 

^n^^-^ = exp(-r^ m,a;/4) , (20) 

where uu is the oscillation frequency, as mentioned earlier. 

Since the Higgs is produced from an S'-wave cc, the perturbative distribution 
amplitude is ultraviolet stable and can be normalized to one, in order to correspond 
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to Pc as a probability to have such a charm quark pair in the proton 

*g(r) = ^i^oKr). (21) 

Here the modified Bessel function KQ{mcr) is the Fourier transform of the energy 
denominator associated with the cc fluctuation. We assume the c and c quarks to 
carry equal fractional momenta. For fixed ag the energy denominator governs the 
probability of the fluctuation in momentum space, since perturbatively one treats the 
charm quarks as free particles. 

Now we can calculate the part of the matrix element in Eq. (13) related to Higgs 
production from the IC. We assume the initial cc wave function to be a linear combina- 
tion (19) of the nonperturbative, Eq. (20) and perturbative, Eq. (21) wave functions, 
and the flnal state wave function of the cc pair in the Higgs in the form (15). The 
result reads, 

j <fTH\T) ^^-'1^ (l - e-^^-') *,,(r) = ^ X^X^ 
V2mc ^ ^ \lmj 



X 



(22) 



— ♦ 

Here we have made use of Mu ^ rric and expanded the exponentials exp(±iA; • r) 
up to the flrst nonvanishing term. We also dropped the integration over a, assuming 
that the c and c in the IC component carry the same momentum, i.e. '^ccij'i «) = 
^cc(?^)^(« — 1/2)- Correspondingly, we have fixed e = Mnjl and we will assume 
1. 

The result of integration in Eq. (22) shows that the perturbative contribution is 
quite enhanced relative to the nonperturbative term. First of all, the enhancement 



by factor ■sj2mc/oJ is due to the fact that the projection to the point-hke Higgs wave 
function is proportional to ^'^(O), and the perturbative fiuctuation has a smaller 
radius. Another enhancement factor, ln(Mj//2mc), is due to the long power tail in 
the momentum distribution in the perturbative IC wave function, which the non- 
perturbative one has a Gaussian cut off. Thus, the perturbative term in the matrix 
element Eq. (22) is relatively enhanced by one order of magnitude. 

Notice that the nonrelativistic nonperturbative solution should not be used for 
convolution with the highly perturbative Higgs wave function. The large transverse 
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momentum tail of the IC should be represented by the perturbative term in Eq. (19). 
Unfortunately, the normalization of such a perturbative tail is unknown, and we 
normalize it to the IC weight 1%. 

Enhancement of the perturbative intrinsic heavy flavor in Higgs production is 
especially large for the top component. Using the IT wave function Eq. (17) we get 



/ 



X 



1 + 



1-5 



ln(l - 5) 



(23) 



where 5 = M^j/Arn^. 

The proton is produced in a similar way from the P-wave 3q in the projectile IC 
component. The exponentials, however, should not be expanded, since the radius is 
not small. Therefore, using the relation, 

d(t) p (e'P-^'^ - e-'P-^'""^ = 47ri p ^ Ji {kp/2) , (24) 



we get for the integral over (Pp in Eq. (13), 



/ 



V^Rsq P 



'■k/2 



TiRr, 



4 Rp Rsq 



e-y [Jo{y) - My)] , 



(25) 



where TZ'^ = 2 R^ Rlg/iR^ + Rjg), and y = k'^/32. For further estimates wc assume 
that Rp — R^q, so TZ — Rp. Since we assumed a meson- type quark-diquark structure 
for the proton, the mean separation R^ — 2(r^^)p/3. The transition proton form 
factor, exp(— A;^7?.^/32), cuts off the integration over (Pk. 

— * 

Now we are in a position to perform the last integration over k in Eq. (13), 

da'^ipp^ppH) 32 GFPjc{z)mt [aZ{s')]' l\s') 



X 



1 + 



dX2 



16 {rl). lis') 
where -f{s') ^ Rl/4Rl{s'). 



1-X2 Mfj Bis') {rl,)p [2 + 27(5') + 7^(5')]' 

2 r I / \ i2 



2m, 



+ Vl - /?2 In 



p 

Mh 
2mr 



(26) 
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The cross section of Higgs production from the intrinsic bottom has the same 
form as Eq. (26), and we assume that the weight of intrinsic heavy flavor scales as, 
PiQ = Picml/mQ. However, as we found above, if the Higgs mass is restricted by 
Mfj < Amf, production from intrinsic perturbative top component of the proton has 
cross section, 



X 



1 + 



dX2 

16 (r,\). 7(s') 



1 — X2 

1-6 



MhB{s') {rl,), [2 + 27(s')+7^(^')P 

-I 2 



1 + 



ln(l - S) 



(27) 



4.2 Feynman Xp distribution of Higgs particles 



The Xp distribution of the cross section Eqs. (26)- (27) is related to the LC wave 
function "ifjciR, z) of the system 3q — cc, namely to the function Pic{z) defined in 
Eq. (8). The momentum fraction z is related to xi^2 and Xp by Eq. (14). The shape 
of Pic{z) strongly correlates with the origin of IC, a nonperturbative component of 
the proton wave function or a perturbative fluctuation. 



4.2.1 Nonperturbative IC 

In principle one can construct hadronic LC wave function by diagonalizing the LC 
Hamiltonian. Here we will use the method of Ref. [43] for the Lorentz boost of the 
wave function, which is supposed to be known in the hadron rest frame. The Lorentz 

boost generates higher particle number quantum fluctuations which arc missed by 
this procedure; however this method works well in known cases [44, 45], and even 
provides a nice cancelation of large terms violating the Landau- Yang theorem [46] . 

We assume that the rest frame IC wave function has the oscillatory form (in 
momentum space), 

4,cW,.) = V^) (^)"*exp . (28) 

Here u stands for the oscillator frequency and fi = MccM^q/ [Mcc+M^q) is the reduced 
mass of the cc and 3q clusters. For further estimates we use Mcc — 3GeV and 
M^q — 1 GeV, although the latter could be heavier, since it is the P-wave. 
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To express the 3- vector Q by the effective mass of the system, Mg// = vQ^ + ^lc+ 



(52 + M|g, one can switch to the LC variables, Q and z, 

Q2 



z{l - z) 



z 1 — z 

Then the longitudinal component Ql in the exponent in (28) reads. 



(Ml - Ml^f Ml + Ml 



3q 



(29) 



(30) 



and the LC wave function acquires the form. 



^ic{Q,z)=K^Pic{z) exp 



8a;// 



{Ml - Ml^Y 



(31) 



where 



X 



Q^(2z-1) 

^2(1-^)2 " 



3/2 



exp 



Ml 



+ 



Ml, 



1 - 



(M|c - M. 



3q 



[i-zy 



(32) 



Now we can calculate the z-dependence of the function Pic{z) defined in Eq. (8), 
which controls the Xi dependence of the cross section, 
Pic{z) 1 da^'^{pp ppH) 1 



(33) 



Pic a^^{pp^ppH) dxi Pic 

This function is plotted in Fig. 3. The distribution sharply picks at z ^ 0.75, as one 
could expect, since the IC pair is heavy and should carry the main fraction of the 
proton momentum. Note that at high energies, in particular at LHC, the momentum 
fraction z coincides with the Feynman xh of the Higgs particle, with a high accuracy 



4.2.2 Perturbative intrinsic heavy flavors 

The light-cone wave function of a perturbative fluctuation p — > \3qQQ) in momentum 
representation is controlled by the energy denominator, 

z(l - z) 



^iq{Q,z,k) oc 



Q^ + z^ml + M^^{l-z) 
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(34) 



6.0 




Figure 3: The distribution of produced Higgs particles over the fraction of the proton 
beam momentum. The dotted, dashed and sohd curves correspond to Higgs produc- 
tion from nonperturbative IC {P — 1), perturbative IC (/3 = 0) and IT, respectively. 

Momentum Q was defined in Fig. 2 and Eq. (5). The effective mass of the QQ 
depends on the intrinsic transverse momentum of the QQ pair, — 4(«;^ + mg). 
It is controlled by the convolution of the IC QQ wave function with the P-wave QQ 
wave function in the Higgs and the one-gluon exchange amplitude (see Fig. 2), which 
has the form. 



dK^ ^iq{Q, z, k) Hqq{k + k/2) - Hqq{k - k/2) 



In 



oc z{l — z) 



Q2+4m^(l-z)+m202 



M|,(l-z)+g2+^2^2 



(35) 



This expression peaks at 1 — 2; ~ mp/Mn, therefore the logarithmic factor hardly 
varies as function of which is restricted by the proton form factor. Making use of 
this we perform integration in Eq. (33) and arrive at the following ^-distribution. 



Nz{l - z) 



In 



\M'ii-iml\{l-z) 
irriQ (1—2) 2^ 



- + 



(36) 



19 



where A^" is a constant normalizing to one the integral over z. The corresponding 
^-distributions for charm and top are shown in Fig. 3 by dotted (dashed) and solid 
curves respectively. 



4.3 Energy dependence 

One can integrate in Eqs. (26)- (27) over xi using relation (14). Since the momentum 
distribution of Higgs produced from the nonperturbative IC sharply peaks at z = 
Zq = 0.75, one can replace Pic{z) =^ 6{z — Zq) Pjc- With a reasonable accuracy we 
can fix z at the same value for the perturbative case and heavier flavors too, which is 
justified by the rather mild dependence on s of other factors in Eqs. (26)-(27). 

Since at high energies 2; ~ ~ 1 — Xi, performing integration in Eq. (26) one 
arrives at, 



a 



IC 



. ^. 32 GfPic m 



X 



1 



zo Mj, B{s) (rl), [2 + 27(5) + ^s)? 



y^+y(r3^in 



2mr 



(37) 



where s = s Zq Mq / Mfj. Analogous expression should be valid for Higgs production 
from intrinsic bottom. For top quark in the proton we use Eq. (27) which leads to. 



a {pp — > ppH) 



tot 



X 



1 + 



16 (r,\). 7(5) 



8 GpPiT m 

Mf, B{s) (rl,), [2 + 27(5) + 72(5)]3 

1 + ^ ln(l - S) 



(38) 



Notice that function 7(5) increases with energy as s^^^^, and such a steep rise of 
the denominator in Eq. (37) is not compensated by the rise of the total cross section 
in the numerator. Therefore, the diffractive cross sections, Eqs. (37)-(38), turns out 
to decrease at asymptotic energies approximately as inverse energy. This unexpected 
result may be interpreted as follows. The source of the falling energy dependence is 
the steep rise with energy of the mean transverse momentum of gluons as is given 
by the unintegrated gluon density Eq. (11), (A;^) = A/Rl{x) oc (s/Af^)°-^*. Also the 
integral over k'^ of the distribution (11) rises with energy, and its value at A; = 
is steeply falling. The rise comes for large transverse momenta which, however, are 
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cut off by tfie nucleon form factor Eq. (20). This is why the diffractive cross section 
(37) is steeply faUing. Indeed, without this form factor, for instance in the reaction 
pp HXp, the cross section would rise as (s/M//)° ''. 

Nevertheless, at the energy of LHC, ^/s = 14 TeV, the effective energy is rather 
low, ^/es = 120 GcV (we assume Mh = 100 GeV) the cross section still rises with 
energy. Indeed, = 0.36 fm^, so 7(es) = 0.55 is still rather small at this energy, and 
the cross sections Eqs. (37)- (38) rise as, 

a'^ipp ^ ppH)lhc a ij^j . (39) 

However, at much higher energies the energy dependence will switch to a steeply 
falling one. Besides, absorptive or unitarity corrections are known to slow down the 
rise of the cross sections. 



4.4 Absorptive corrections 

The amplitude of any off-diagonal large rapidity gap process is subject to unitarity 
or absorptive corrections, which have the intuitive meaning of a survival probability 
of the participating hadrons. To include these corrections one should replace the 
diffractive amplitude as, 

s) ^ /f|(6, .) [1 - Im fJib, s)] , (40) 

The data for elastic pp scattering show that the partial amplitude /^f (&, s) is 
constant energy at small impact parameters fe — 0, while rising as function of energy 
at large b [47, 48, 49] . This is usually interpreted as a manifestation of saturation of 
the unitarity limit, Im/^^ < 1. Indeed, this condition imposes a tight restriction at 
small 6, where Im/^^ ^ 1, leaving almost no room for further rise. We will treat the 
Pomeron as a Regge pole without unitarity corrections: 



,2 



2B:fis) 



(41) 



where aZ{s) = 21.8 mb x {s/M^y, and e = 0.08; B^J{s) = B% + 2a'jp ln(s/M2) with 
— 7.5GeV~^. Due to the accidental closeness of 2a']p/B'li — 0.067 and e, the 
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pre-exponential factor in (41) hardly changes with energy even without unitarity cor- 
rections. It is demonstrated in Ref. [49] that not only at 6 = 0, but in the whole range 
of impact parameters, the model Eq. (41) describes correctly the energy dependence 
of the partial amplitude /gf (&, s). 

Thus we arrive at the absorption corrected cross section, 



u^Q(j)p ppH) = a^^{pp — > ppH) 



X 



1-1 



+ 



TT B{s') + 25J(.0 (47r)2 Bl^is') [B{s') + i?J(.')] 



(42) 



This is not a severe suppression even at the energy of LHC, where the absorptive 
factor is 0.2. 

Including the absorptive corrections we calculated the total cross sections for 
diffractive Higgs production, pp Hpp, from the IQ components. The results at 
the energy of LHC, ^/s = 14TeV, are plotted as function of Higgs mass in the 
Fig. 4. We assume a pcrturbative origin for all intrinsic components, a 1/itLq scaling 
for their weights, and 1% probability of IC ior P — in Eq. (37). Note that the 

1 



10 



a 10 

(X 

g 1 



10 
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Figure 4: The cross section of the reaction pp — > Hp + p as function of the Higgs 
mass. Contributions of IC (dashed) of IB (dotted) and IT (solid). 



contributions of the intrinsic charm and bottom falls steeply with the mass of the 
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Higgs in accordance with Eq. (37). The contribution of the intrinsic top rises with 
Mh unless Mh > 2mt ~ 350 GeV, then the cross section starts falhng. 

5 Further possibilities to get a larger cross section 

5.1 Direct production of Higgs from a colorless IQ 

A heavy flavor QQ pair in the IQ component of the proton may be found in a colorless 
state. In this case the Higgs particle can be produced directly from this pair via 
Pomeron exchange as is shown in Fig. 5. We consider the example of intrinsic charm 



H(r) 




Figure 5: Higgs production via Pomeron exchange 

of a nonperturbative origin throughout this section. 

At first glance one may think that this channel is less suppressed by powers of the 
Higgs mass compared to the mechanism presented in Fig. 2. Indeed, the integration 
over k does not have the upper cutoff imposed by the proton form factor in the 
previous case; therefore it may compensate two powers of Mh in the amplitude. This 
analysis is correct. Nevertheless the amplitude turns out to be more suppressed than 
in the diagram Fig. 2. 

The diffractive amplitude A{x2,Pi,P2) is proportional to the matrix element of 
the dipole cross section crqq{r) between the initial cc wave function "^ccir) and the 
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distribution amplitude of the cc in the Higgs, 

A(x2,pi,P2) oc J d\H\r)a-,g{r)^-^{r) . (43) 

This factor contains all the dependence on the Higgs mass. To estimate it one can 

use a Gaussian shape for both ^ccl'") and H{f). Then one finds A oc ^/m^/Mjj. A 
more refined calculation confirms this, 



2mcUJ Rq{s) \2' ' ArricUJ 



(44) 



where the initial state cc pair is assumed to be in a P-wave. Here U {a, b, x) is the 
confluent hypergeometric function, and we use its asymptotic behavior at x ^ 1, 

U{a, b, x) = x-"- + Oix-"-^) . (45) 

Notice that a convolution similar to Eq. (43) also defines the scale dependence of 
the amplitude of photoproduction of heavy quarkonia, which also behaves like (Q^ + 

Af^)^'^. Thus, the complementary mechanism of diffractive Higgs production, besides 
the smaller probability to find the colorless IC component, is additionally suppressed 
by 1/Mfj compared to the dominant mechanism depicted in Fig. 2. Therefore, this 
contribution can be safely neglected. 



5.2 Nuclear enhancement 

The produced Higgs is supposed to escape detection and to be identified only by using 
the missing mass spectrum. One may also consider the same reaction on a bound 
proton in pA collisions where the nuclear debris spectators flying in the same direction 
as the Higgs. The nuclear enhancement in this case is not as large as one could naively 
expect. The reason is that absorptive corrections are stronger than those considered 
above in Sect. 4.4 for the case of pp coUisions. The survival probability represented 
by the last factor in Eq. (42) can be evaluated within the Glauber approximation for 
pA collisions as 

Zeff^jJd%T4b)e-^^T.W, (46) 
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where Ta(&) = dz pA{b, z) is the nuclear thickness function at impact parameter 

— * 

6, and PaQ>) is the nuclear density. We assume here that diffractive recoil neutrons 
cannot be detected; otherwise the factor (46) should be multiplied by A/Z. 

The nuclear enhancement for lead according to (46) is rather mild, of about an 
order of magnitude, since Ze// ~ 2.5 should be compared with the suppression factor 
0.2, for pp collisions calculated in Sect. 4.4. Gribov corrections [50] are known to 
make nuclei more transparent, therefore they may substantially increase the survival 
probability factor [51, 31]. If we employ the simplest quadratic dependence of the 
dipole cross section cr^q(r) oc r^, then the nuclear enhancement is considerably larger, 
a factor of about 50. 

6 Conclusions 

The key assumption underlying our analysis of high xp Higgs hadroproduction is the 
presence of intrinsic heavy flavor \uudQQ > Fock components in the proton bound- 
state wave function. Such quantum fluctuations are in fact rigorous consequences of 
QCD. The the probability for intrinsic heavy quark on the heavy quark mass falls 
as -j§r- in non-abelian theories and can be computed from the operator product 
expansion [11]. In such Fock states the heavy quarks Q and Q carry the highest 
light-cone momentum fractions. Thus although they have small probability, intrinsic 
heavy quark Fock states are highly efficient in transferring the momentum of the 
proton to the momenta of particles in the final-state, especially heavy quarkonium 
and the Higgs which can sum the momenta of both the Q and Q. It is thus interesting 
and important that measurements of the production of heavy quarkonium at high xp 
as well as other heavy hadrons such as the and Ac be carried out at RHIC, the 
Tevatron, as well as the LHC in order to test this novel feature of QCD. 

As we have reviewed in section 3, there is substantial but not conclusive phe- 
nomenological evidence for intrinsic charm at the 1% probability level in the proton. 
It is thus particularly important to have measurements of the charm and bottom 
structure functions in deep inelastic lepton-proton scattering over the full range of 
Xfej. One must allow for intrinsic sea components at any scale Qq when parameter- 
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izing the proton's structure functions, since the intrinsic Fock states are responsible 
for the u{x) ^ d{x), s{x) ^ s{x) asymmetries, as well as the high-x c{x) and b{x) 
distributions. There are also important nuclear and heavy quark threshold effects 
related to intrinsic charm which can be tested at lower energy fixed-target facilities 
such as JLAB, GSI-FAIR, and J-PARC [52]. 

As we have emphasized here, the materialization of intrinsic heavy flavor states 
in the proton leads to Higgs production in the proton fragmentation region: this 
includes inclusive production pp — > HX, singly diffractive production pp — > p + HX 
and exclusive diffractive production pp — > p + H + p, reactions which should be 
considered in addition to the conventional central rapidity production processes. The 
fractional momentum distribution for a Higgs produced by combining the momenta of 
both heavy quarks in the IQ Fock states is presented in Fig. 3. As seen in the figure, 
the Higgs can be produced with momentum fractions as large as xp ~ 0.8 or even 
higher. One also produce the Higgs inclusively from leading-twist PQCD processes 
such as gc — > He and gb — > Hb where the high momentum of one intrinsic heavy 
quark is transferred to the Higgs. 

We have focused in this paper on diffractive exclusive Higgs production pp — > 
p + H+p, since in principle, only the final-state protons need to be measured and the 
Higgs can be reconstructed from the missing mass distribution. We note, however, 
that detecting the diffractive signal pp ^ p + H + p poses new challenges: When the 
Higgs is produced at large xp, one of the final-state protons will be need to be detected 
at a small momentum fraction 1 — xp, which is outside of the usual acceptance of 
forward proton detectors. 

The underlying color structure of the intrinsic Fock state and the gauge theory 
properties of the two-gluon exchange mechanism for high energy diffraction play key 
roles in the physics of the exclusive diffractive Higgs hadroproduction process. The 
main result of our analysis, the cross sections given in Eqs. (37,38), demonstrates that 
the heavier the intrinsic heavy quark, the larger is the cross section for the doubly 
diffractive reaction. It rises with itlq linearly if the heavy quarks are confined by a 
potential, and is presented in Fig. 4 if the QQ appear in the proton as a perturbative 
fluctuation. The production cross section also steeply rises with energy oc s'^^'^, which 
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is characteristic for the energy dependence of hard reactions. Absorptive correction 
slow down this rise and eventually stop it at very high energies, above the energy 
range of LHC. Asymptotically, in the Proissart regime, this cross section is expected 
to fall. Numerical predictions for diffractive Higgs production from IC, IB and IT 
components arc shown in Fig. 4. The cross section will be further enhanced from 
possible Fock states of the proton containing supersymmetric partners of quarks or 
gluons. We also have discussed a potential increase in the rate for such reactions 
using proton-nucleus collisions. 
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